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SUMMARY 

By the use of an integral technique , the laminar boundary- layer 
equations are reduced on the wind-ward generator of the plane of symmetry 
to a set of simultaneous algebraic equations. The Chapman-Rube sin 
* temperature- viscosity relation and a Prandtl number of 1 are assumed. 

The method enables the skin friction coefficients and Stanton number to 
be calculated in a much shorter time than was needed to obtain exact 
numerical solutions from the boundary- layer equations. The solutions ob- 
tained by this method are, for the most part, within 5 percentage points 
of the exact solutions. 
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INTRODUCTION 

Some designs for supersonic aircraft and missiles indicate the use of 
a conical nose or forebody. For design purposes it Is necessary to know 
the skin friction and heat transfer on the conical surface for all pos- 
sible flight conditions. One important condition Is flight at angle of 
attack, since the boundary- layer flow is complicated by the addition of a 
crossflow, which increases with Increasing angle of attack. 

Several analytical papers have been published dealing with the effect 
of angle of attack on the boundary layer of a cone. Laminar flow over 
cones at small angle of attack has been studied (ref. 1) by using a lin- 
earization of the nonlinear boundary-layer equations for a cone as derived 
in reference 2. In references 3 and 4, the results of reference 1 have 
been extended to cover the combined effect of spin and small angle of 
attack. In reference 5, the boundary- layer equations for large angles of 
attack were solved by a numerical method but the solutions were li m ited 
to the plane of symmetry. The preceding solutions were limited to the 
insulated case with a Prandtl number equal to 1. 
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Additional numerical solutions were computed in reference 6 for the 
equations for large angle of attack in the plane of symmetry. Heat- 
transfer solutions were also obtained in reference 6 by noting that for 
a Prandtl number of 1 the energy equation in the plane of symmetry could 
be written in the same form as the momentum equation along the generator 
and, therefore, the enthalpy profile is identical to the velocity 
profile .in this direction. An approximate correction factor for heat 
transfer is given in reference 6 for the cas_e of a Prandtl number not 
equal to 1. 

Obtaining a numerical, solution to the boundary- layer equations of 
reference 5 is an involved process and requires an electronic computor. 
It was therefore desirable to find an approximate method of solving the 
boundary- layer equations from which results could be obtained with only 
the use of a desk calculator. An integral method is presented herein, 
that results in a pair of simultaneous algebraic equations in place of 
the nonlinear ordinary differential equations of reference 5 for the 
plane of symmetry. 

Interpolation of the exact solutions of reference 6 gives good re- 
sults within the range of variables studied^ but the method developed 
herein enables quick calculation of results outside the range of varia- 
bles considered in reference 6. For this purpose, a computational pro- 
cedure for the present method Is given in appendix B. The technique 
given in this paper also indicates a method for studying the boundary 
layer on the cone other than In the plane of symmetry. 


EQUATIONS OF MOTION 

The boundary-layer equations are given in reference 2 for a general 
orthogonal coordinate system in which the body surface is defined by 
y = 0 and the coordinate x is defined so that the body cross sections 
are similar for various values of x. These equations can therefore be 
applied to a cone at angle of attack with the coordinate x along a 
generator, y perpendicular to the cone surface, and xf3<p on the cir- 
cumference of the cone, where <P is the meridional angle measured from 
the windward generator on the plane of symmetry and (3 is the sine of 
the cone semivertex angle £ . Figure 1 shows this coordinate system. 
Since the flow past a cone at angle of attack is still conical, there is 
no pressure gradient along a generator of the cone. Therefore, for 
steady- state conditions the boundary- layer equations from reference 2 
can be written as 
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(lb) 


(le) 


(l d) 

(le) 

Equations (la) and (lb) are momentum equations in the x and q> 
directions, respectively, equation (lc) is the continuity equation, equa- 
tion (id) is the energy equation, and equation (le) is the equation of 
state. All symbols are defined in appendix A. 

The folio-wing physical assumptions, -which are basically the same as 
in reference 5, are used in this paper: 

(1) A thin laminar boundary layer across which the static pressure 
is constant and also Cp is constant 

(2) Frandtl number equal to 1 

(3) The Chapman- Rubesin temperature- viscosity relation (ref. 7) is 


dw , dv 

pU 'Sx + pV By 


bw pw dw puw _ 1. bp b f dw\ 

3y xp x “ " xp &p dy v 'by) 

1 d(pw) 


pu 


^ 

xp 


d(pu) + 5(pv) 

dx ' ay xp d<P ‘ x 
dE , dE . pw dE /du , dv , 1 dw u\ 

-s + PT -s^- , -Sp^p=-i , ('s + ^ + iF^ + 5) 

r— _ . 




JL 

^e 



( 2 ) 


where the constant C is evaluated to match the Sutherland value of vis- 
cosity at the cone surface 


C = 



I 1 / 2 T e + 198.6 
' T g + 198.6 


(3) 


(4) The cone surface temperature is constant. This differs from the 
assumption of zero heat transfer given in reference 5 in that heat trans- 
fer is considered. 


For Pr = 1, equations (l) can be combined to give the new energy equation 
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/ 5h ^ Sh , w Sh\ 5 / 5h\ 

P \ U- 5x v "Sy xp "3<p/ "Sy V 1 Sf) 

where 

H = c^T + ^ (u 2 + w 2 ) 

If, in addition, the surface temperature is constant and a new 
quantity ® is defined as 

H - EL 
0 = ^ 


the energy equation (4) can be written as 

The boundary conditions on the cone at angle of attack are 
u=v=w-0=O at y=0 

and 


/ S0 ^ as ^ 
p \ u 3SZ + V 3y + 


w as\ a ( a®\ 


(*) 


(5) 


( 6 ) 


(7) 


(8) 


u = U e , V = V e = 0, w = W e , @ = 1 at y = & 

where y = S at the outer edge of the boundary layer. It is assumed 
that both the velocity and thermal boundary layers have the same thick- 
ness 5. 


Combining each of the momentum and energy equations with the con- 
tinuity equation and then integrating across the boundary layer from 
y = 0 to y = 6 results in 

^ uu + _L ^uw + ^uw d(pgu|) ^ £uu 

ax xp a cp Bxo u 2 a cp x 

H K e e 


1 

xBu e dcp 



dy + - 
x 



■a *2 
p w 
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P~e 

p e u e 



wJ 


y=0 


(9a) 
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d£-WTl + _1_ + £wV _ d(p e uf) ^ ^WU 

Sx X8 d9 x8p e u| X 


1 

xP u e 


3"e f * * 1 

-3J / P w dy-- 
•'0 

^@u 1 

~ 5x~ xp ^<p 


pW dy 


^Qw 5 (Pe u e) 
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8 dp 
' *Bp e u| 

+ 

x " Pe u e 



where 


t = 
b uu 

/•B 

/ p*u*(l - u*)dy 

'o 

(10a) 

t, 

b uw ^ 

/•B 

/ p*v*(l - u*)dy 
0 

(10b) 

^ 

/* 5 

/ p*v*(v* - w*)dy 

b 

(10c) 

^wu = ^ 

f p*u*(v| - w^Jdy 

b 

(lOd) 

^0u = j 

/*8 

f p*u*(l - 0*)dy 

0 

(lOe) 

^0w = J 

/»5 

f p*v*(l - 0*)dy 

(lOf) 


The starred quantities introduced into equations (9) and (10) are dimen- 
sionless quantities formed as follows: 

u* = u/u e p* = p/p e 

w* = w/u e (i* = ^/h e 

T* = T/T e 


©* = 0 


(H) 
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Equations (9) and (10) are reduced to incompressible form, by means of the 
transformation 



( 12 ) 


The relation between the boundary- layer thickness in the physical 
coordinates 8 and the thickness in the transformed coordinates A was 
obtained from equations (5), (6), and (12) . The pressure gradient in the 
# direction was evaluated at the outer edge of the boundary layer. By 
using the Chapman-Rubesin temperature-viscosity relation (2) and the 
dimensionless variable 





the transformed boundary-layer equations become 

Se uu 1 ^uv . 9 uw d(p e ui) . 9 UU 


A ^e 
xpu e "5W 



drj 



dT] 


Pe^ ( 5u*\ 

p e u^ v^r 



(13) 


(14a) 


(14b) 
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d%u 1 <^@w 00V d(p e Ue) e 6u _ PeC f d®*\ 

~W + Xp ~W + Xpp e ii e d# X = P e u e^ \^T /t^q 


where 



and 


1 _ u e 


Ti " 2c p T e 

The boundary conditions are 

u* = V* = 0* = 0 at T] 

and 


* _ * n * *1 
U =0 = 1, V = v p 


du dw 


d0 


at TJ a* 1 


dr] dt] dT] 


= 0 


(14c) 

(15a) 

(l^b) 

(15c) 

(l5d) 

(l5e) 

(I5f) 

(16) 

(17) 
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APPROXIMATE VELOCITY AND ENTHALPY PROFILES 

Approximate solutions for equations (14) using the boundary con- 
ditions given by equation (17) can be obtained by assuming a polynomial 
to represent the boundary-layer profile for each of the velocity and 
thermal boundary layers. The parameters in_the assumed profile are eval- 
uated from the boundary conditions. Both a third- and fourth- degree 
polynomial were assumed. The results from the fourth-degree polynomial 
differed more from the exact solutions of reference 6 than those using 
the third-degree polynomial. Since it is desired to use an approximation 
that agrees as closely as possible to the exact solutions in the plane of 
symmetry , the method using the third- degree polynomial Is presented. 

Assume that 

u* = a^T} + &2 r l 2 + a 3 1 l 3 
w* = b^q + bgT] 2 + b 3 t) 3 

0 = C-j^T] + CgTj + CgT) 

Substituting boundary conditions (17) into equations (18), the 


assumed profiles become 

u* = (3 - 2 t))t} 2 + aii(l - 2i) + tj 2 ) (19a) 

w* = (3 - 211 ) 1 ^* + bi](l - 2q + t i 2 ) (19b) 

©* = (3 - 2ii)t) 2 + cq(l - 2q + T] 2 ) (19c) 

The expressions given in equations (14) and the integrals given in 

equations (15) become 

°uu - § [| + ! « - 1 * 2 ] < 2 °*> 

- § [l + ! k - 1 b2 ] < 2ob ) 


(18a) 

(18b) 

(18c) 
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J 0 ^ dT l = 35 L r ~ * 3 . 

] (20J) 

uV di] . i [13 TfJ + awj + if b + 

i abj (20k) 

(ir^o = 8 

(20Z) 

(|£) = 1 

7r P o 

(20m) 

(iH-o ’ ° 

(20n) 
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The "boundary- layer thickness is assumed to develop according to the 
expression 

A " xl/2 .. (21) 

where K is an unknown proportionality parameter. The quantities K, 
u*, p_, and u are assumed to vary with # hut not with X. The 

assumption of expression (21) for the laminar "boundary- layer thickness 
is based on a statement in reference 5 that Blasius-type parabolic sim- 
ilarity exists in meridional planes. 

Substitution of equations (20) and (21) into equations (14) gives 
three ordinary differential equations in four unknowns: a, b, c, and 

c/k 2 . 


KANE OF SYMMETRY 


The equations of motion for the boundary layer can be greatly sim- 
plified by restricting them to the plane of symmetry. In the plane of 
symmetry # is equal to either 0 or it, however, both theory and ex- 
periment indicate that the boundary-layer assumptions do not hold for 
4> = at large angles of attack, since the boundary layer either becomes 
too thick or separation occurs. Therefore, the boundary- layer equations 
will be solved herein only for the case $ = 0 where, because of sym- 
metry, w = w e = w* = 0 and therefore b = Oj also p'(<$) = 0. 

Inspection of the ^-momentum equation (14b) shows that each term 
vanishes on the plane of symmetry. In order to use this equation, it is 
necessary to differentiate it with respect to $ before restricting it to 
the plane of symmetry. For w = 0, equations (14a) and (14c) are iden- 
tical so that u* = ©* and a = e. It is, therefore, only necessary 

to solve the X-momentum equation and the 4? derivative of the 
$ -momentum equation. 


The resulting equations in the plane of symmetry are 





35a ~ (22a) 
K 2 
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(15 - x ) l 1 + ^)(f - § a )] ■ - 3S Jr f < 2Zb > 

These equations are simultaneous algebraic equations in the three un- 
knowns: db/d<£, C/K 2 , and a. 


Examination of the exact boundary- layer profiles as given in ref- 
erence 6 indicated similarity in the x-direction that was practically 
independent of angle of attack, temperature ratio, and cone angle. There- 
fore, a is assumed constant for all values of the variables »nd its value 
can be determined by comparison with the exact solutions. 


The quantities 



are functions of the flow at the outer 


edge of the boundary layer and can be evaluated from data tabulated in the 
M.I.T. tables of flow over a cone (refs. 8 and 9). As given in reference 
5, these quantities are 



where cl is the angle of attack, 8 is the sine of the cone s emi vertex 
angle, and all other quantities on the right side of the equations are in 
the notation of references 8 and 9. The barred quantities refer to zero 
angle of attack. 


The quantities given in equations (23) and (24) are exact up to the 
order of a 2 but terms of higher order of a are considered negligible. 
This assumption sets an upper limit on the angle of attack. This limit 
is also a function of the other variables of the problem. 
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COEFFICIENTS OF FRICTION AND HEAT TRANSFER 


The expressions for the components of viscous shear stress for the 
direction along the most windward generator and the circumferential di- 
rection are 



\ p u 2 
2 e e 



< j >=0 




0 



(25a) 


(25b) 

(25c) 


In terms of the quantities obtained from this report , equations (25) can 
be -written as 


( C f,x\>=0 “ 2a K 1 ^ 2 

(26a) 

( c f,<p v=° = 0 

(26b) 

2 — — (Re )- V 2 

\H = o ■ K a# 

(26c) 

Similarly, the expression for the Stanton number is 


St ' <>eM H a S ‘ 


X s (f Lo 

" PeV H as ' V 

(27) 


which, in terms of the presented. Quantities becomes 

St = a | (Re x ) _1 / 2 


( 28 ) 
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SOLUTIONS OF EQUATIONS 

The exact solutions for the boundary-layer equations as given in 
reference 6 were used to evaluate the constant a. The expressions given 
in reference 6 for ski n -friction coefficients and Stanton number are 


( c f>x}p = o 2f w 

(29a) 

_ 3 p kl(r rt 73 C 

\ <^<P 4=0 w V Re x 

(29b) 

(£ H=o = Jr^ 

(29c) 


where f" , T|r", and ® * are defined in reference 6 and 
w w w 

k = -!('£:£) 

- 3(5 Vdi 4 =0 


(30) 


It should be noted here that in reference 6, as in the present paper , the 
energy equation was similar to one of the momentum equations and therefore 

■ ®w- 

The skin- friction coefficients as given in equations (26) were set 
equal to the exact skin- friction coefficients as given In equations (29) 
for the particular case of k = 0, which is for zero angle of attack. 
Combining these coefficients with the X-momentum equation., the resulting 
value for a Is 


a = 1.5921 

Using this value of a, the boundary- layer equations for the wind- 
ward generator (eqs. (22a) and (22b)) become 


„ _ . 2.7753 dw e . 1.3026 db _ _ C 
7.2737 + — — aT + — — a* = 55.724 ^ 


(31a) 
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di> 


|1 ( dw eY l dv i 


db 1 

d<f» P 


8.5000 



+ 1.4167 


dw£ 


e db 

d4> d# 


+ 0.6667 

CL'S* 


db 


P \d# 


Y 1 \( l\ 

/ + Vd Wj 35 + 17 * 706 t7 + 12.856 ^ 


C db 

- - 35 £2 25 


(31b) 


These boundary-layer equations can be solved simultaneously to give 
db/d# and C/K? once the conditions of the problem are given. 


Solutions of the boundary- layer equations for the windward generator 
are given in table I along with the exact solutions of reference 6 for _ 
those conditions for which exact solutions have been obtained. The solu- 
tions are given in terms of the quantities obtained in reference 6. The 
relations between these quantities and the quantities presented herein are 



( c f,x\p=0 

2 



0.9192 



(32a) 


1 Y^f/p\ .px 0.5774 db (IT 
^w“ 3pk \ / Q 7 % «f\ d$ \ k 2 

\“**A-o 



(32b) 


(32c) 


The data presented in table I are plotted in figure 2. 

The errors between the approximate solutions obtained by the present 
method and the exact numerical solutions of reference 6 are less than 1 
percent for both Stanton number and skin- friction coefficient along the 
most windward generator. With the exception of the extreme cooling case^ 
where I b /Tq a* 0; the errors are less than 5 percent for the derivative of 
the circumferential skin- friction coefficient. In the range of conditions 
covered in reference 6, interpolation of the _exact solutions gives better 
results than the use of the present method. However, in regions not 
covered by the exact solutions the present method will give useable results. 
A computational procedure is given in appendix B. _ 
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The integral method used reduces the boundary- layer equations to 
ordinary differential equations over those parts of the cone where the 
boundary-layer assumptions hold. Figure 2 shows that the solutions of 
these equations, when reduced to the plane of symmetry, agree closely 
with exact solutions of the original equations. Therefore, solutions of 
the ordinary differential equations away from the plane of symmetry can 
be obtained and will give some idea of the boundary-layer flow off the 
plane of symmetry. 


CONCLUSIONS 

By means of an integral technique, the nonlinear partial differential 
equations for the laminar boundary layer over a cone at angle of attack 
are reduced in the plane of symmetry to a set of simultaneous algebraic 
equations. Third-degree polynomials are assumed for both velocity and 
enthalpy profiles. The approximate solutions agree with exact solu- 
tions for a Prandtl number of 1 within 1 percent for Stanton number and 
skin friction coefficient along the most windward generator and, for the 
most part, within 5 percent for the derivative of the circumferential 
skin-friction coefficient. 

The method presented gives an easy method of determining skin- 
friction coefficients and Stanton numbers outside the range for which 
exact solutions are available. 

Due to the close agreement between the approximate solutions and the 
exact solutions on the plane of symmetry, the same technique can be used 
to study the boundary- layer flow off the plane of symmetry. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 4, 1958 
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APPENDIX A 
SYMBOLS 


a l* a 2 jS 3 

b 

Tdi,Td 2> 13 3 

C 

C f 

c 


V 

c l’ c 2 } °3 
E 


f" 

v 


H 

K 


k 

Pr 

P 

q. 

R 

Re x 


parameter In generator velocity profile 

coefficients in assumed generator velocity profile 

parameter in circumferential velocity profile 

coefficients in assumed circumferential velocity profile 

constant in Chapman-Rube sin temperature-viscosity relation 

skin-friction coefficient 

parameter in enthalpy profile 

specific heat at constant pressure 

specific heat at constant volume 

coefficients in assumed enthalpy profile 

Internal energy, c v T 

generator shear parameter (defined in ref. 6) 
enthalpy as defined in equation (5) 

proportionality parameter in transformed boundary-layer 
thickness "~ 

circumferential velocity- gradient parameter 

Prandtl number 

static pressure 

heat- transfer coefficient 

gas constant 

Reynolds number based on distance x 


St 


Stanton number 
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u 

V 

"W 


X 


X 

Y 

y 

a 


P 

A 

5 


n 

0 


i 

e iJ 

8 ' 

v 


p- 

p 

$ 

<p 


temperature 

dimensionless external temperature 

velocity in x- direction 

velocity in y-direction 

velocity in circumferential direction 

transformed coordinate along generator 

coordinate along generator 

transformed coordinate perpendicular to surface of cone 

coordinate perpendicular to surface of cone 

angle of attack 

sine of cone semivertex angle 

boundary-layer thickness in transformed coordinates 
boundary-layer thickness 
defined by equation (10) 

dimensionless coordinate perpendicular to cone 
enthalpy difference ratio 
semi vert ex angle of cone 
defined by equation (15) 

heat- transfer parameter (defined in ref. 6) 
coefficient of thermal conductivity 
viscosity 
density 

transformed angular coordinate 

angular coordinate measured from most "windward generator 
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i|f^. circumferential shear- gradient parameter (defined in ref. 6) 

Subscripts : 

as adiabatic surface 

e outer edge of boundary layer 

s surface 

x at distance x along generator 

* 

cp at angular distance tp 

0 free- stream stagnation 
Superscripts : 

* dimensionless quantity 

1 differentiation with respect to independent variable 


'4854 
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APPENDIX B 


COMPUTATIONAL PROCEDURE 


The analysis given in the present paper enables one to calculate 
skin- friction coefficients and Stanton numbers on the windward generator 
of a cone at angle of attack. The a n alysis is, however, limited to the 
case of Prandtl number equal to 1. It is necessary to know the following 
conditions : 


(1) Mach number upstream of shock wave 

(2) Cone surface to free-stream stagnation temperature ratio 

(3) Semi vertex angle of cone 

(4) Angle of attack 

The following quantities are to be evaluated: 



where a and 8 = sin A are given and the other quantities on the right 
side of the equations are given in references 8 and 9. The barred quan- 
tities refer to zero angle of attack. 


Substitute the values of (dWg/d$^_Q and 1 /t^ into the following 
equations : 

„ _ , 2.7753 dw t . 1.3026 db _ _ C 
7.2737 + — = 55.724 
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These equations can be solved simultaneously to obtain db/d# and 
C/Y ? . Since the second equation is quadratic in db/d#; two solutions 
will occur. The desired solution is that which makes C/K? positive. 
The skin- friction coefficients and Stanton numbers are then obtained 
directly from 
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TABLE I. - COMPARISON OF APPROXIMATE AND EXACT SKIN- 
FRICTION AND HEAT- TRANSFER PARAMETERS FOR MOST 
WINDWARD GENERATOR OF A CONE AT ANGLE OF 
ATTACK FOR FRANDTL NUMBER OF 1 



i 

k 

or 0^ 

* 

tt 

w 

Exact 

Approx- 

imate 

Exact 

Approx- 

imate 

0 

0 

0 

0.3321 

0.3321 

0.4238 

0.4514 



.6 

.433 0 

.4318 

.5527 

.5837 




.5143 

.5124 

.6570 

.6923 


2.5 

0 

0.3321 

0.3321 

0.6532 

0.7464 



.6 

.4598 

.4563 

.8596 

.9673 



BC 1 ' 

.5569 

.5516 

1.0281 

1.1518 


5.0 

0 

0.3321 

0.3321 

0.8826 

1.0407 



.6 

.4815 

.4764 

1.1422 

1.3249 




.5898 

.5823 

1.3634 

1.5739 

m 

0 

0 

0.3321 

0.3321 

0.5923 

0.6079 



.6 

.4468 

.4455 
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Circumferential velocity-gradient parameter, k 


(a) Ratio of surface temperature to free-stream total temperature, T s /Tq ;= 0. 

Figure 2. - Shear and heat-transfer parameters from exact and approximate solutions in plane of symmetry. 
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Circumferential velocity- gradient parameter, fc 


(h) Ratio of surface temperature to free- stream total temperature, T s /T q n 0.5. 

Figure 2. - Continued . Shear anil heat- transfer parameters from exact ami approximate solutions in plane 
of symmetry. 
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Circumferential velocity-gradient par a meter, i 

(c) Ratio of surface temperature to free- stream total temperature, T a /Tg m 1.0. 

Figure 2. - Concluded . Shear and heat- transfer parameters from exact and approximate solutions In plane 
of symmetry. . 
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